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Abstract 

The dynamics of a spatially extended system of two competing species in the presence of two 
noise sources is studied. A correlated dichotomous noise acts on the interaction parameter and a 
multiplicative white noise affects directly the dynamics of the two species. To describe the spatial 
distribution of the species we use a model based on Lotka-Volterra (LV) equations. By writing 
them in a mean field form, the corresponding moment equations for the species concentrations 
are obtained in Gaussian approximation. In this formalism the system dynamics is analyzed for 
different values of the multiplicative noise intensity. Finally by comparing these results with those 
obtained by direct simulations of the time discrete version of LV equations, that is coupled map 
lattice (CML) model, we conclude that the anticorrelated oscillations of the species densities are 
strictly related to non-overlapping spatial patterns. 
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I. INTRODUCTION 



The dynamics of real ecosystems is strongly affected by the presence of noise sources, such 
as the random variability of temperature, resources and in general environment, with which 
the system has a multiplicative interaction . In this paper we analyze the time evolution 
of a spatially extended system formed by two competing species in the presence of two noise 
sources. We get the dynamics in the formalism of the moments. We study the role of the two 
noise sources on the ecosystem dynamics, described by generalized Lotka-Volterra equations 
in the presence of external fluctuations, modelled as multiplicative noise. Specifically we 
focus on the time behavior of the I s * and 2 nd order moments of the species concentrations. 
We find that the 1 st order moments are independent on the multiplicative noise intensity. 
On the other hand the behavior of the 2 nd order moments is strongly affected by the presence 
of a source of external noise. We find anticorrelated time behavior of the species densities. 
Comparing our results with those obtained by calculating the same quantities within a 
coupled map lattice (CML) model |3j, we conclude that the anticorrelated oscillations of 
the species concentrations are strictly related to non-overlapping spatial patterns 14]. Our 
theoretical results could match data from a real ecosystem, whose dynamics is affected by 
the random variability of the^environment, and could provide useful tools to predict behavior 
of biological species [l|, 




II. THE MODEL 

Our system is described by a time evolution model of Lotka-Volterra equations, within 
the Ito scheme, with diffusive terms in a spatial lattice with iV sites 
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m,j = mA 1 - w - P x i,j) + yi,jV°ytlj + D $^(2/7 - ViJ), ( 2 ) 
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where x it j and y^ denote respectively the densities of species x and species y in the lattice 
site is the growth rate, D is the diffusion constant, and indicates the sum over 

all the sites. Here ^fj(t) and £,fj{t) are statistically independent Gaussian white noises with 
zero mean and unit variance, a x and a y are the intensities of the multiplicative noise which 
models the interaction between the species and the environment, and j3 is the interaction 
parameter. 
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A. The interaction parameter 

Depending on the value of the interaction parameter, coexistence or exclusion regimes 
take place. Namely for (3 < 1 both species survives, while for f3 > 1 one of the two species 



extinguishes after a certain time. These 



Volterra's deterministic model 



DOB 



;wo regimes correspond to stable states of the Lotka- 



]. Moreover periodical and random driving forces 
connected with environmental and climatic variables, such as the temperature, modify the 
dynamics of the ecosystem, affecting both directly the species densities and the interaction 
parameter. This causes the system dynamics to change between coexistence {(3 < 1) and 
exclusion (/3 > 1) regimes. To describe this dynamical behavior we consider as interaction 
parameter (3{t) a dichotomous stochastic process, whose jump rate is a periodic function 

7(0 

*)={\ , n (3) 
I 7q (1 + A | cosu;t|) , At > r d . 

Here At is the time interval between two consecutive switches, and is the delay between 
two jumps, that is the time interval after a switch, before another jump can occur. In eq. Q, 
A and u = (2tt)/T are respectively the amplitude and the angular frequency of the periodic 
term, and 70 is the jump rate in the absence of periodic term. This causes (3{t) to jump 
between two values, Pdown < 1 and (3 up > 1, which correspond to the dynamical regimes of 
the deterministic Lotka- Volterra's model (coexistence and exclusion regions). Because the 
dynamics of the species strongly depends on the value of the interaction parameter, we report 
in Fig.[T]the time series of f3{t) for different values of delay r^, namely = 10, 43.5, 60, with 
Pdown = 0.94 and f3 up = 1.04. We note that the correlation time of the dichotomous noise 
affects the switch time between the two levels of (3{t). For a delay time a bit less than T/2, 
we observe a synchronization between the jumps and the periodicity of the rate 7(t). This 
synchronization phenomenon is due to the choice of the value, which stabilizes the jumps 
in such a way they happen for high values of the jump rate, that is for values around the 
maximum of the function 7(t). This causes a quasi-periodical time behavior of the species 
concentrations x and y, which can be considered as a signature of the stochastic resonance 
phenomenon in population dynamics 0, Q, 1^ • Therefore we fix the delay at the value 
td = 43.5, corresponding to a competition regime with /3 switching quasi-periodically from 
coexistence to exclusion regions (see Fig. [TJd). 
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FIG. 1: Time evolution of the interaction parameter (3(t) with initial value /3(0) = 1.04 and 
different values of delay: (a) = 10 (a), 43.5 (b), and 80 (c). The values of the other parameters 
are: A = 9.0, cj/(2tt) = 10~ 2 , 70 = 2 • 10~ 2 . 



III. MEAN FIELD MODEL 

In this section we derive the moment equations for our system. Assuming N — > 00, we 
write Eqs. ([TJ and (j2J in a mean field form 

x = f x (x, y) + ^g x (x)C + D(< x > -x), (4) 
y = f v (x, y) + ^g y {y)i v + D(< y >- y ), (5) 

where < x > and < y > are average values on the spatial lattice considered (strictly 
speaking they are the ensemble average in the thermodynamics limit) and we set f x (x, y) = 
fix(l -x-Py), g x (x) = x, f y (x, y)=fiy(l-y- fix), g y (y) = y. By site averaging Eqs. (gj) 
and ((SJ), we obtain 

<x> = <f x (x,y)>, (6) 
<V> = <f y (x,y)>. (7) 

By expanding the functions f x (x,y), g x (x), f y (x,y), g y (y) around the 1 st order moments 
< x > and < y >, we get an infinite set of simultaneous ordinary differential equations for 
all the moments [10j . To truncate this set we apply a Gaussian approximation, for which 
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the cumulants above the 2 order vanish. Therefore we obtain 



< x > = /i<x> (1- < x > -(3 < y >) - fi(/3fiu + a*2o) (8) 

< y > = /i < y > (1- < y > -(3 < x >) - fi((3fi n + M02) (9) 
/i 20 = 2/i/i 20 - 2D/i 20 - 2/i(3 < y > fi 2 o ~ 2\i < x > (/3/in + 2fi 20 ) 

+ 2a x {< x > 2 +/i 20 ) (10) 
(iQ2 = 2/i/i 02 - 2Dfi 02 - 2///3 < x > /i 2 - 2/i < y > (f3fx n + 2fi 02 ) 

+ 2ct v (< j/ > 2 +/102) (11) 
/in = 2/i/in - 2D fi n - <x> [2fjLfjLn + (jlP(jj,u + ^02)] 

- < y > [2/x/in + /i/3(//u + ^20)] , (12) 

where /i 2 o, A*o2> A*n are the 2 nd order central moments defined on the lattice 

H20 = < x 2 > - < x > 2 (13) 

H02 = <y 2 > -<y> 2 (14) 

fin = < xy > — < x >< y > . (15) 



In order to get the dynamics of the two species we analyze the time evolution of the I s * and 
2 nd order moments according to Eqs. (jffi-(|12|). We fix the delay time at the value Td = 43.5, 
corresponding to a quasi-periodic switching between the coexistence and exclusion regimes, 
and we obtain the time series of the moments for two values of the multiplicative noise 
intensity a = o x = a y , namely a = 10 -12 , 10 -6 , and in the absence of it. The values of 
the parameters are // = 2, D = 0.05. The initial values of the moments are < x(0) > 
= < 2/(0) > = 0.1, /i2o(0) = ^02(0) = /in(0) = 0. These initial conditions correspond 
to uniformly distributed species on the lattice considered. In Fig. El we note that the 1 st 
order moments of both species oscillate together quasi regularly around 0.5, independently 
on the multiplicative noise intensity (see Fig. |2k). The noise intensity affects strongly the 
dynamics of the 2 nd order moments. In the absence of noise H20, 1^02, A*n are zero. For very 
low levels of multiplicative noise (er = 10~ 12 ) quasi-periodical oscillations appear with the 
same frequency of the interaction parameter f3(t), because the noise breaks the symmetry of 
the dynamical behavior of the 2 nd order moments (see Fig. |2b) Q • The time behavior of the 
variances of x and y species coincides all the time with alternating periods, characterized 
by small (close to zero) and large values. However the negative values of the correlation fin 
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FIG. 2: Time evolution of the I s * and 2 nd order moments. The time series of (a) < x(t) > 
and < y(t) >, and (b) fi2,o and /io,2 respectively, are completely overlapped. The values of the 
multiplicative noise intensity are: a = 0, 10~ 12 , 10~ 6 , from top to bottom. Here = 43.5, fi = 2, 
and D = 0.05. The initial values of the moments are < x(0) > = < y(0) > = 0.1, fi2o(fy = 
^02(0) = /Un(0) = 0. The values of the other parameters are the same of Fig. ^ 



indicate that the two species distributions are anti-correlated. This means that the spatial 
distribution in the lattice will be characterized by zones with a maximum of concentration 
of species x and a minimum of concentration of species y and viceversa. The two species 
will be distributed therefore in non-overlapping spatial patterns. This physical picture is 
in agreement with previous results obtained with a different model jq|. For higher levels 
of multiplicative noise (a = 10~ 6 ) the amplitude of the oscillations increases both in /j, 2 o, 
fi02 and fin. This gives information on the probability density of both species, whose width 
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and mean value undergo the same oscillating behavior. The anti-correlated behavior is 
enhanced by increasing the noise intensity value (see Figs. |2t>) • We note that the amplitude 
of the oscillations in Fig. |2b increases with the noise intensity a and it is of the same 
order of magnitude. The periodicity of these noise-induced oscillations shown in Fig. |2] is 
the same of the interaction parameter (3(t) (see Fig. [TJ. Even if it is due to a very different 
mechanism, this behavior is similar to the stochastic resonance effect produced in population 
dynamics, when the interactionparameter is subjected to an oscillating bistable potential in 
the presence of additive noise [3,0]. We note that in the absence of external noise (er = 0) 
both populations coexist and the species densities oscillate in phase around their stationary 
value 7]. This occurs identically in each site of the spatial lattice. The behavior of the mean 
value therefore will reproduce this situation. For a ^ 0, anticorrelated oscillations appear 
due to the multiplicative noise, superimposed to the average behavior obtained for a = 
and distributed randomly in the spatial structure. By site averaging these noise-induced 
oscillations (see Ref. p|) we recover the average behavior obtained in the absence of noise. 
This explains why the first moment behavior is independent on the external noise intensity. 



IV. COUPLED MAP LATTICE MODEL 

In order to check our results we consider a different approach to analyze the dynamics 
of our spatial extended system. We consider the time evolution of CML model, which is 
the discrete version of the Lotka-Volterra equations with diffusive terms. For this model we 
found anticorrelated spatial patterns of the two competing species p|, that are related to the 
dynamical behavior of the moments of the species densities. Here we calculate the moments 
in the CML model. Within this formalism, the dynamics of the spatial distribution of the 
two species is given by the following equations 

= " " ^M?) + y/wflg!* + D £>W - (16) 

7 

*/S +1) = ^S } (! - y% - (B) *S } ) + v^/S^f + D - vth ( 17 ) 

7 

where x\j and ylj denote respectively the densities of prey x and prey y in the site 
at the time step n, /x is the growth rate and D is the diffusion constant. and £]j are 
independent Gaussian white noise sources with zero mean and unit variance. The interaction 
parameter (3^ corresponds to the value of (3{t) taken at the time step n, according to Eq. Q. 
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Here ^ 7 indicates the sum over the four nearest neighbors. To evaluate the I s * and 2 order 
moments we define on the lattice, at the time step n, the mean values < x >( n > , < y >^ n ' , 

' 1 ' CML a CML 



( n ) 



the variances vari , var y 



uarW = Vs* =V — - — , z = x } y (19) 

and the correlation coefficient corr^ n > of the two species 



c°rr <"> = (20) 



(«) 



with 

co <) = — - ^ : : (2i) 



E^(^S ) -<^> (w) ) (yS ) - < v > (n) ) 



The number of lattice sites is N = 100 x 100. The time behavior of these quantities, for two 
levels of the multiplicative noise and in the absence of it, is reported in Fig. El The 1 st and 
2 nd order moments of Eqs. ()18j) - f!21l) correspond respectively to the same quantities shown 
in Fig. |21 We note that the two set of time series are in a good qualitative agreement. The 
discrepancies in the oscillation intensities are due to the different values of the stationary 
values of the species densities in the considered models. Specifically: x st = y s t = a/(l + /3) — 
0.5 for the mean field model, and x n st — y™ t = (1 — l//x)/(l + (3 n ) ~ 0.25 for the CML model. 
Moreover the behavior of the 2 nd order moments in Fig. Eb shows little irregularities with 
respect to that obtained in the mean field model, because the species interaction in the CML 
model is restricted to the nearest neighbors. 



V. CONCLUSIONS 

We report a study on the dynamics of a spatially extended ecosystem of two competing 
species, described by generalized Lotka-Volterra equations. Two noise sources are present: a 
multiplicative white noise, which affects directly the two species densities, and a correlated 
dichotomous noise, which produces a random interaction parameter whose values jump 
between two levels. The role of the dichotomous correlated noise is to control the dynamical 
regime of the ecosystem (see Fig. while the multiplicative noise is responsible for the 
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FIG. 3: (a) Mean values: < x > CML , < y > C mli an d (b) variances: var x , var y and corr 
of the two species, as a function of time. The values of the multiplicative noise intensity are: 



a = 0, 10 ,10 , from top to bottom. The initial values of the species concentrations are x 
Vij = 0-1 f° r au s hes The values of the other parameters are the same of Fig- El 



(0) 



anticorrelated behavior of the species concentrations (see time behavior of /in in Figs. EJd 
and Eh). The anti-correlated oscillations are enhanced by increasing the multiplicative noise 
intensity. The mean field approach with the Gaussian approximation enables us to obtain 
the time behavior of the I s * and 2 nd order moments, which characterize the spatio-temporal 
behavior of the ecosystem. We compare the results obtained within a mean field approach 
with those obtained with a CML model. The agreement is quite good and allows us to 
conclude that the spatial patterns of the two species, within the mean field approach, should 
be non-overlapping as those obtained with the CML model . Our theoretical results could 
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explain the time evolution of populations in real ecosystems whose dynamics is strictl" 
dependent on noise sources, which are always present in the natural environment 0, 
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